Coexistence of amplitude and frequency modulations in intracellular calcium dynamics 
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The complex dynamics of intracellular calcium regulates cellular responses to information encoded 
in extracellular signals. Here, we study the encoding of these external signals in the context of the 
Li-Rinzel model. We show that by control of biophysical parameters the information can be encoded 
in amplitude modulation, frequency modulation or mixed (AM and FM) modulation. We briefly 
discuss the possible implications of this new role of information encoding for astrocytes. 

PACS numbers: Valid PACS appear here 



Many cells use calcium signaling to carry informa- 
tion from the extracellular side of the plasma membrane 
to targets in their interior This information serves 
many different purposes, from triggering the developmen- 
tal program of fertilized mammalian eggs to mediating 
neural activity and learning or inducing cell death. Since 
virtually all cells employ a network of biochemical reac- 
tions for Ca 2+ signalling, much effort has been devoted 
to understand the functional role of the Ca 2+ response 
and to decipher how this complex dynamical response is 
regulated by the biochemical network of signal transduc- 
tion pathways @, [§]. About a decade ago, several experi- 
ments indicated that Ca 2+ signals in response to external 
stimulus can encode information via frequency modula- 
tion (FM) or in some other cases via amplitude mod- 
ulation (AM) Consequently, it has been shown that 
these observations can be captured separately by minimal 
models consisting of two dynamical variables such as the 
Li-Rinzel [f| or the Dupont-Goldbeter models @. It was 
also shown that higher-order models (ones with several 
dynamical variables and/or intracellular diffusion mecha- 
nisms) can exhibit different and more advanced encoding 
modes 0. 

Here we propose that under certain conditions, hetero- 
geneous dynamics of intracellular Ca 2+ could also be ex- 
plained by opportune parameter modulation within min- 
imal models. More specifically, we employ arguments 
of bifurcation theory to illustrate that within the mini- 
mal Li-Rinzel model the same cell could encode the in- 
formation about external stimuli in amplitude modula- 
tion (AM) of calcium oscillations, in frequency modu- 
lation (FM) or in both (AFM). Our work is motivated 
by the calcium signalling in astrocytes, a predominant 
non-neuronal (glial) cell type that plays a crucial role 
in the regulation of neuronal activity [g, 0] . We explain 
why for this case, our results can be crucial for a bet- 
ter understanding of synaptic information transfer and 
propose that they might be equally important for better 
understanding of other examples of processes regulated 
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TABLE I: Parameters used in the original Li-Rinzel model. 

by Ca 2+ signalling. 

Calcium dynamics is controlled by the interplay of 
calcium-induced calcium release, a nonlinear amplifica- 
tion process regulated by the calcium-dependent open- 
ing of channels to Ca 2+ stores such as the endoplasmic 
reticulum (ER) , and by the action of active transporters 
(SERCA pumps) which enable a reverse flux. The dy- 
namical variables of LR model, that is studied here, are 
the free cytosolic Ca 2+ concentration (C), and the frac- 
tion of open inositol trisphosphate (IP3) receptor sub- 
units, h: 
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The dynamics of C is controlled by three fluxes, corre- 
sponding to: 1. a passive leak of Ca 2+ from the ER to 
the cytosol, (Ji ea k)] 2. an active uptake of Ca 2+ into ER, 
Jpump, due to the action of the pumps; 3. a Ca 2+ release 
{J chan) that is mutually gated by Ca 2+ and the inositol 
trisphosphate (IP3) concentration, (J): 



Jleak(C) = V 2 (C -(l+C 1 )C) 



K\ + C 2 
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(3b) 



Jchan(C,I) = v 1 m 3 oa h 3 (C -(l + c 1 )C) (3c) 
The gating variables and their time-scales are given by: 
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with Q 2 = d 2 ■ The level of IP 3 is directly controlled 
by signals impinging on the cell from its external environ- 
ment. In turn, the level of IP3 determines the dynamical 
behavior of the above model. One can therefore think 
of the Ca 2+ signal as being an encoding of information 
regarding the level of IP 3 . 

The original set of biophysical parameters, as given 
in table HI corresponds to AM encoding. For these pa- 
rameters, the phase-plane and bifurcation analysis re- 
veals that, at I « 0.355/iAf, limit-cycle Ca 2+ oscillations 
emerge through a supercritical Hopf bifurcation. From 
fig-Ilk h is evident that the amplitude of Ca 2+ oscillations 
increases between the two bifurcation points, while fig. [lb 
shows that the frequency of the oscillations is almost con- 
stant - hence the term "Amplitude Modulation" . Am- 
plitude modulation by IP3 has been observed in many 
experiments; however new findings indicate that, under 
some conditions, variations (by external stimulation) in 
the level of IP3 can also lead to frequency modulation 
[loj ]. These observations motivated us to re-examine the 
LR model to investigate if changes in the biophysical pa- 
rameters could lead to frequency-modulated dynamics. 
A nonlinear system can exhibit frequency modulation in 
the presence of saddle- node bifurcation [ill ] . This lat- 
ter describes a transition of a system from an excitable 
state (in which there are three fixed points: stable, un- 
stable and a saddle) to a limit cycle. At a certain value 
of the control parameter / = I sn , the stable and saddle 
fixed points coalesce and the only remaining attractor is 
a limit cycle. The frequency of the oscillations of the 
limit cycle can be very sensitive to the distance from the 
bifurcation point (i.e.: I — I sn )i whereas the amplitude 



remains almost constant [12 1 



We have explored the range of parameters for which 
the LR system can exhibit a saddle-node bifurcation with 
the level of IP3 being the control parameter. We found 
that K3 (the affinity of the active SERC A pump) , d$ (the 
receptor affinity for IP 3 ) and v 2 (the rate of Ca 2+ leakage 
from the ER) all can regulate the switching between AM 
and FM encoding dynamics. We further discovered the 
existence of a new dynamical regime in which the vari- 
ations of the IP3 are co-encoded both in amplitude and 
frequency modulations ("AF Modulations"). This AFM 
dynamics exists for higher levels of the cell-averaged rest- 
ing Ca 2+ concentration Co as well as for lower U3 rates 
of Ca 2+ uptake by SERC As. We present a biophysical 
picture of these different regimes and comment on the 
physiological implications of these results with particular 
attention to astrocytes. Although there have been some 
earlier indications that the LR model can encompass ex- 
citable behavior ljj 14 1, these works did not present a 
complete analysis nor a biophysical picture. For brevity 
we consider only the case of varying K3. 

We begin with the well studied and simpler AM 
dynamics that corresponds to higher values of K3. As 
stated above, in this case a limit cycle emerges through a 
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FIG. 1: The Li-Rinzel model, (a) Bifurcation diagram for 
the original set of parameters of the Li-Rinzel model: (— ) sta- 
ble fixed points, (• • • ) unstable ones, (•) stable limit cycles, 
(o) unstable ones. Oscillations are born via supercritical Hopf 
bifurcation at [IPs]~ 0.355 /xM and die via subcritical Hopf 
bifurcation at [IPs]~ 0.637 /iM. While the amplitude changes, 
the frequency is nearly constant (b). (c) Nullclines (green: h, 
orange: C) for the case of an unstable point, (d) At basal IP3 
levels (~ 0.015 fiM.) J pU m P (red curve) intersects J re ; ( — ) at 
a calcium-level such that J' rel < 0. This situation also occurs 
at higher IP3 when J re i becomes bell-shaped. 
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FIG. 2: An excitable version of the Li-Rinzel model. 

(a) Bifurcation diagram and (b) period diagram of an ex- 
citable version of the LR model with K3 = 0.051 /itM. In 
this case, four bifurcations exist: a saddle-node and a saddle- 
node on invariant circle (SNIC) at [IPs]~ 0.479 fj,M and 
[IP3]~ 0.526 j«M, and two subcritical Hopf bifurcations at 
[IP 3 ]~ 0.51 (UM and [IP 3 ]~ 0.857 /M. (c) Between the two 
saddle-node bifurcations nullclines intersect in three points 
which are a stable focus (•) and an unstable node (□) sepa- 
rated by a saddle (A), (d) J pU m P (magenta curve) now inter- 
sects Jrei at lower Ca 2+ . 



supercritical Hopf bifurcation where a single stable fixed 
point becomes unstable. In fig. [TJ: we show the nullclines 
of h and C when the fixed point is unstable. In this case 
the properties of the limit cycle can simply be understood 
from linear stability analysis of the unstable fixed point 
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near the bifurcation. In fig. QJl we plot the calcium fluxes 
(for two different values of 7) as determined by setting 
h to /100(C): these fluxes capture the fast time scale re- 
sponse of the system close to the fixed point, since the 
rate of receptor response is slower than that of the Ca 2+ 
concentration. The fixed point occurs at high calcium, 
where the slope of the efflux J re i — J c han + Jleak is neg- 
ative. At values of Ca 2+ above the unstable fixed point, 
Jrei < J pump, and there is a return flow in the (h,C) 
phase plane towards small values of C and vice versa be- 
low the fixed point. The instability, then, gives rise to a 
dynamical flow that does not exhibit any strong positive 
amplification and instead is due to the slow dynamics. 
This fixes the frequency to be close to that of the h time 
delay which does not vary much across the oscillatory 
regime. Conversely, the amplitude is relatively free to 
vary so that the system exhibits amplitude modulation 
in response to IP 3 variations. 

From the above analysis, it is clear that the key to get 
frequency modulation is to make the stable fixed point 
occur at low calcium, on the rising part of the efflux 
curve. In our case this is accomplished by increasing the 
pumping rate by considering lower values of K3 , as shown 
in fig. [2] In such conditions in fact, when we consider 
the h and C nullclines (fig. [2fc) we note that for proper 
I values there is a stable fixed point at low Ca 2+ levels 
which is close to a saddle-point. Inspection of the char- 
acteristics of the Ca 2+ fluxes (shown in fig. [2ji) reveals 
that at the saddle-point the slope of the characteristic 
of J re i is steeper than that of J 
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In this situation, a 
finite yet relatively small deviation away from the stable 
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FIG. 3: AFM encoding (colors online), (a) Continuation of 
AM and FM bifurcations allows to identify a region comprised 
between a Bautin (B) and a cusp (CP) bifurcation where 
AFM-encoding could be found ((□) Hopf points, (A) saddle- 
node points), (b) Modulation map illustrating the regime 
of existence of the three classes - AM, FM and AFM - of 
dynamical responses, (c) Bifurcation diagram and (d) pe- 
riod diagram for an AFM-encoding version of the LR model 
(d 5 = 0.2 fM, Co = 4.0 fM). 



fixed point crosses the saddle-point separatrix and leads 
to a large excursion in the phase plane. For this reason, 
the Hopf bifurcation of the stable fixed point (which still 
occurs before I sn for the new set of parameters) must 
now be subcritical and different deviations (i.e. different 
values of I) lead to trajectories with similar amplitudes 
(fig.[2^L), as this is determined by the global flow. At the 
same time the flow field in the vicinity of the stable fixed 
point and the saddle-point is very weak, hence the period 
of the excursions is very sensitive and can be effectively 
modulated by the IP3 (fig. This dynamics shows 

"Frequency Modulation" . 

The transition from AM to FM occurs via a charac- 
teristic codimension-2 bifurcation sequence which com- 
prises the following steps: first, the lower supercritical 
Hopf point changes into a subcritical one via Bautin bi- 
furcation; then, elsewhere in the parameter space, a cusp 
bifurcation generates the SNIC and saddle-node bifur- 
cations which are responsible for variable-period oscil- 
lations (fig. EK)- When the Hopf bifurcation is not yet 
strongly subcritical but we can nonetheless feel the influ- 
ence of saddle-node coalescence, we can predict that the 
emerging oscillations would show significant variability 
both in amplitude and frequency (figs. [3l:,d, Hfc). This is 
the previously mentioned AFM dynamics where IP3 vari- 
ations modulate both the amplitude and the frequency 
of the oscillations. 

According to our analysis, AFM encoding can be found 
in several cases, most typically for higher Co (fig. [3Jd) 
or smaller V3 values (table III)) . In terms of the Ca 2+ 
fluxes (eqs. [3]), we note that such a choice of parameters 
counteracts the effect of a lower A3 by increasing the 
distance between the characteristics, either by increasing 
J re i (through an increase of Co) or by reducing J pU mp 
(by lowering V3). It follows that the stable fixed point 
moves towards higher calcium levels and the amplitude 
modulation can coexist with frequency modulation. 

In summary, we showed the existence of three distinct 
classes of information encoding modes in the response 
of Ca 2+ to IP 3 variations - AM, FM and AFM. Two 
of these classes, AM and FM, were previously demon- 
strated in other minimal models whereas AFM has 
been explained so far only by extended models which 
included either diffusive terms 0] or complex feedback 
loops @, H, [TBI. We found that AFM dynamics could 
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TABLE II: Parameter ranges for the coexistence of amplitude and 
frequency modulation of Ca 2+ response in Li-Rinzel model. 
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FIG. 4: Different types of excitability. Proper tuning of 
parameters allows for different Ca 2+ responses for a generic 
IP3 stimulus (depicted in (d)). (a) The original LR parame- 
ters (table[I} provide amplitude variability of oscillations that 
occur at almost fixed frequency, (b) A higher SERCA pump 
Ca 2+ affinity (K3 = 0.051 fiM) is responsible for oscillations 
with variable frequency but nearly constant amplitude. De- 
spite the different nature of bifurcations underlying these two 
cases, both AM and FM features can also be found for exam- 
ple in presence of higher cell-averaged total free Ca 2+ levels 
as shown in (c) for Co = 4 fiM (with K 3 = 0.051 ^M). Stim- 
ulus: (a) AIP3 = 0.4, 0.1, 0.1, 0.5 fjM; (b) AIP 3 = 0.4, 0.2, 
0.2, 0.5 fiM; (c) AIP 3 = 0.17, 0.03, 0.09, 0.5 /jM. 



also be reproduced by a minimal model for some partic- 
ular range of the model parameters. These findings hint 
that by activating intracellular mechanisms that control 
the values of physiological parameters that correspond to 
the model parameters C n or v 2 , K 3 and v 3 for SERCA 
pumps or ds for IP 3 Rs [16j, the type of information en- 
coding can be regulated. 

We note that for the different modes of Ca 2+ response 
to have an information encoding role, a corresponding 
decoding mechanism must exist. The existence of decod- 
ing mechanisms of AM and FM have been proposed in 
171 ] based on model studies of the cooperative binding 
of Ca 2+ to generic effector enzymes. In principle, the 
same mechanism can also decode information embedded 
in Ca 2+ that corresponds to AFM. 

In the context of communication theory, cellular Ca 2+ 
signalling can be regarded as a bifurcation-based encod- 
ing system: the baseline IP3 level Io is set to be suf- 
ficiently close to a bifurcation point so that the varia- 
tions in IP3 caused by external signals regularly cross 
that point. This mechanism can be in principle realized 
in new kind of electronic systems. In AM, Ca 2+ peaks 
encode the information; in FM, variations in the IP3 will 
trigger bursts of Ca 2+ spikes (fig. [4} with information 
encoded in the inter-spike intervals. In the mixed AFM 
mode, both features carry information which can be sep- 
arately decoded by different downstream effectors with 
different Ca 2+ responses. This is particularly suitable 
for those systems that require special constraints in co- 



ordination of informational input from multiple channels, 
for which recently, mixed AM and FM modulation (MM) 
has been receiving much attention. 

The above can be very relevant for the case of astro- 
cytes regulation of synaptic information transfer. As- 
trocytes respond to synaptic activity through their in- 
tracellular Ca 2+ dynamics which in turn feeds back to 
neurons by triggering release of gliotransmitter. AFM 
encoding could have deep consequences [3]. We might 
expect that the short-time scale effectors (mostly sensi- 
tive to the number of pulses) are involved in feedback to 
the local synapse whereas the long-time scale ones (which 
integrate the total signal) coordinate information with 
other astrocytes via intercellular signalling. 
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